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We introduce a strong-disorder renormalization group (RG) approach suitable for investigating 
tiie quasiparticle excitations of disordered superconductors in which the quasiparticle spin is not 
conserved. We analyze one-dimensional models with this RG and with elementary transfer matrix 
methods. We find that such models with broken spin rotation invariance generically lie in one 
of two topologically distinct localized phases. Close enough to the critical point separating the 
two phases, the system has a power-law divergent low-energy density of states (with a non-universal 
continuously varying power-law) in either phase, due to quantum Griffiths singularities. This critical 
point belongs to the same infinite-disorder universality class as the one dimensional particle-hole 
symmetric Anderson localization problem, while the Griffiths phases in the vicinity of the transition 
are controlled by lines of strong (but not infinite) disorder fixed points terminating in the critical 
point. 



I. INTRODUCTION 

Recently, there has been considerable theoretical activ- 
ity concerning the effect of disorder on the quasiparticle 
spectrum inj-dirty superconductors with different pairing 
symmetries .Eli3 The basic philosophy underlying most of 
these developments is to start with a weakly disordered 
problem and investigate the fate (in some RG sense) of 
disorder at large length scales using field theoretic meth- 
ods. B An analysis of this sort depends most crucially on 
the possible symmetries of the quasiparticle Hamiltonian, 
and it is the corresponding universal properties that have 
received most attention thus far. These include the lead- 
ing effects of disorder on the conductivity and density of 
states in the delocalized regime, as well as the universal 
scaling properties of the localization transition. 

Our focus in this article is quite different: We ask if 
Grifhths effects, whereby rare configurations of the disor- 
der potential over large regions of space give rise to non- 
universal contributions that dominate some low-energy 
property, may be important in dirty superconductors. In 
particular, are there situations in which such Griffiths ef- 
fects lead to a singular enhancement in the quasiparticle 
density of states near the Fermi energy? Unfortunately, 
this intriguing possibility has not attracted much atten- 
tion in previous work. 

Such issues are difficult to address within the weak- 
disorder framework mentioned above. Instead, we in- 
troduce a RG approach that is suitable for situations in 
which the effective value of disorder becomes large at 
low energies. Our basic result is simply stated: It is 
indeed possible for such Griffiths effects to dominate the 
low-energy properties, at least in one-dimensional models 
with either triplet pairing, or strong spin-orbit effects, or 
spin-flip scattering off frozen magnetic impurities. More 
specifically, we show that such models with broken spin 
rotation (SR) invariance generically lie in one of two 



distinct localized phases. These two phases are distin- 
guished by a topological property (reffected in the pres- 
ence or absence of zero-energy end-states in a large but 
finite wire) that makes it impossible to go smoothly from 
one phase to the other. At the phase transition separat- 
ing the two, the effective disorder grows without bound 
when viewed on ever-smaller energy scales, and the low- 
energy j-density of states behaves as p{e) ~ l/e|ln'^e| 
(DysorO singularity). The localized phases near this 
transition, on both sides, are Griffiths phases having a 
power-law singularity (with a non-universal exponent) in 
the low-energy density of states; these Griffiths phases 
are themselves characterized by strong (but not infinite) 
effective disorder in the limit of low energies. [Note that a 
recent weak-disorder analysis, Refs. || and ^, of this prob- 
lem has been carried out only at criticality. The results 
of Refs. ||j6| for the thermal conductance and the den- 
sity of states are consistent with the predictions of our 
RG approach at such a critical point — however, as we 
show here, the geneac behavior of the system is localized 
rather than critical.Q] 

The detailed scaling properties of the low-energy 
strong-disorder critical point and the nearby Griffiths 
phases are, of course, specific to our one-dimensional ex- 
amples. However, our RG approach is well suited for 
studying possible Griffiths effects in two or more dimen- 
sions as well, and some speculations along these lines are 
briefiy discussed towards the end of this article. 

II. PHYSICAL PICTURE AND MOTIVATION 

Before plunging into the details of our analysis, it is 
useful to have a heuristic picture of the basic physics 
of our Id problems. To this end, we consider a sim- 
ple toy model for a disordered spinless (more physically, 
spin-polarized triplet) superconducting 'wire' with a sin- 
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gle transverse mode (channel) active. We model our sys- 
tem by a lattice Bogoliubov-de Gennes (BdG) Hamilto- 
nian 



H 



Ve„cJ,c„. (1) 



In this toy model, the nearest neighbor hopping ampli- 
tude t and the "p-wave" pairing amplitude A are real 
constants, while the impurity potential e„ can take on 
values Vi and V2 with some probabilities p and 1 — p re- 
spectively. Furthermore, we stipulate that \Vi \ < Vc and 
1^2! >K, where K = 2|i|. 

The significance of the critical value Vc is readily seen 
by solving the pure problem with fixed e„ = V, i.e., with 
the chemical potential equal to —V. One easily finds 
that there is a gap in the quasiparticle excitation spec- 
trum around the Fermi energy for both \V\ > Vc and 
\V\ < Vc- However, at the critical point \V\ = Vc, the 
system is gapless. Thus, there are two different gapped 
phases for our pure system (the phase with \ V\ > Vc sim- 
ply corresponds, in the absence of the pairing term, to 
a situation in which the Fermi level has gone below the 
bottom of the band or above the top of the band; thus, it 
is essentially a 'band insulator'). For our purposes here, 
the important distinction between the two phases has to 
do with low-energy bound states at the ends of a long but 
finite wire of length L with free boundary conditions. In 
the 'gapped superconductor' phase with \ V\ < Vc, such a 
wire has a single quasiparticle state below the gap with an 
excitation energy that is exponentially small in the length 
L; the corresponding wavefunction has weight only in the 
vicinity of the two ends of the chain (in the language of 
Ref. 1^, in the L — > 00 limit, we thus have two zero-energy 
Majorana fermions, one at each end of the chain) . In the 
other gapped phase, \ V\ > Vc, there is no such low-energy 
quasiparticle state.l3 

Now, imagine a disorder realization in which the po- 
tential has value Vi throughout the region between, say, 
sites and L, and value V2 out to a distance L on either 
side of this central segment. The value of the potential is 
left unspecified in the rest of the system. The probability 
of this happening is p^{l — p)^^ = e~^^ , with c appro- 
priately defined. Now, the central region can be thought 
of, for large L, as a finite wire in the phase with \ V\ < Vc 
surrounded by vacuum (this is reasonable since the long 
segments on either side can be roughly thought of as re- 
gions with no particle because the Fermi energy —V2 has 
gone below the bottom of the band; effective couplings 
of the central region with the rest of the system, medi- 
ated through such isolating segments, are exponentially 
small in L). Such a situation will lead to a low-energy 
quasiparticle state with excitation energy el ~ e~''^, 
with some c' of order one. Since such low-energy states 
living on 'domain walls' between large regions in 'oppo- 
site' phases can happen anywhere along the entire length 
of our wire, such disorder configurations will give a non- 
zero contribution in the thermodynamic limit to the den- 
sity of states. This contribution can be estimated as 



PGrisie) = JdL 5{e - e-'='^)e-=-^ - e-i+V^, where we 
have introduced the dynamical exponent z = c'/c. This 
serves as a lower bound on the actual low-energy density 
of states in our disordered problem — thus, we generically 
have a power-law behavior (with a non-universal contin- 
uously varying exponent) of the density of low-energy 
excitations in our toy model. Although this model is 
admittedly crude, the picture of rare configurations of 
disorder over large regions of space leading to singular 
low-energy behavior is at the heart of the more precise 



strong-disorder RG analysis of Section VI 



We conclude this section with some comments on our 
choice of toy-model, and, more generally, on the results 
obtained in this article. Firstly, note that we completely 
ignored the self-consistency condition that, in iwrnciple, 
determines A in terms of the other parameters. 1130 This 
is not expected to matter; in fact, the precise choices 
made for various parameter values are not very impor- 
tant for our conclusion — nor is it important that our toy- 
model has time reversal symmetry. As far as this model 
is concerned, the only important thing is the existence of 
two different gapped phases, with one of them support- 
ing zero-energy end-states. [Note that this is the main 
distinction of our models without SR-invariance from the 
models with SR-invariance: When the quasiparticle spin 
is a good quantum number, there are no such end-states 
in any finite open chain, and we expect no Griffiths ef- 
fects in this case — this will also become clear from our 
more detailed RG analysis.] 

Another concern is that we are treating a one- 
dimensional superconductor with the BdG equations, 
which ignore quantum fluctuations of the condensate 
order parameter A. When the superconductor is in 
more than one dimension, A does have a nonvanishing 
static component at zero temperature, but in a strictly 
one-dimensional system, divergent quantum fluctuations 
mean that the superconducting state does not have true 
long range order or a gap. However, in highly anisotropic 
quasi-one-dimensional superconductors, Tc (and the gap) 
can, in principle, be large compared to the interchain 
hopping energy, and our approach should then apply in 
the range between these two energy scales (while no such 
regime appears to exist in quasi-one-dimensional super- 
conductors known so far, a priori, there are no physical 
reasons that would prevent this from happening in some 
casest3) . Another possible physical realizationEl is that of 
a vortex in a three dimensional gapped superconductor 
in the presence of frozen magnetic impurities or spin- 
orbit scattering. In such a situation, the effect of disor- 
der on the quasiparticle states confined to the vortex core 
can be analyzed by a Id BdG equation approach such as 
the one we employ. Naturally, the choice of probability 
distributions for various quenched random variables will 
be different depending on the physical realization one is 
interested in (for instance, it is more natural to use a 
quenched random A with zero mean when considering 
the vortex problem). However, as will be clear from our 
later analysis, the precise form of the probability distri- 
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bution for various bare couplings in the problem does not 
play an important role in determining the nature of the 
low-energy physics. 

Finally, note that the effects of the residual quasipar- 
ticle interactions are beyond the scope of our analysis. 



III. FORMALISM, SYMMETRY CLASSES, AND 
OUR ID MODELS 

It is^pLiseful to set up notation and review some 
basicsEjcl before proceeding to our actual calculations. 
To this end, consider a general lattice BdG Hamiltonian 



i Im(t+A) 



Q/3 



(2) 



where we use a composite label a = {i,fJ-} for the site 
and spin indices of fermion orbitals. Hopping amplitudes, 
the effects of spin-orbit interaction on the hopping am- 
plitudes, spin- flip scattering from frozen magnetic impu- 
rities, and random potential terms corresponding to non- 
magnetic impurities are now all included in ft-i^jv, while 
pairing amplitudes are represented by A^^^^v. Hermitic- 
ity requires hpa = h*^i3, and we choose A^jq, = — A^^ 
consistent with the fermion anticommutation relations. 
Additional restrictions (to be reviewed below) arise when 
T-invariance is a good symmetry (we will not consider 
cases with SR-invariance in this article). 

The spectrum of quasiparticle excitations for Hsc is 
specified by the spectrum of a (Hermitian) matrix 



n = 



h A 

-A* -h* 



(3) 



I2N (where N is the number of lattice 



acting in an enlarged "particle-hole" Hilbert space. 
A particle/hole mixing unitary transformation C/q = 

\/2 1^1 i 

sites), which acts independently on states correspond- 
ing to each a, transforms Ti into an antisymmetric pure 
imaginary form 

nin.-Uo ^Uo^ ,Im(/i-A) )■ 

This representation is well-suited for a discussion of 
heat transport properties of the quasiparticles. Indeed, 
the quasiparticle thermal conductivity is simply propor- 
tional to fcsTcr, where a is the T — conductivity of a 
(normal) system of non-interacting fermions described by 
the lattice Schroedinger equation corresponding to Tiim 
(see Ref. || and references therein). [Note that Him may 
also be obtained, as in Ref. ||, from the original Hsc by 
writing everything in terms of Majorana fermions and 
then doubling the system, and we will therefore use a 
"copy index" K below to label the different blocks of 

Him-] 
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a) 



b) 



FIG. 1. a) General pure imaginary representation of hop- 
ping (t) and pairing (A) couplings between two fermion or- 
bitals a and /3; b) in the T-invariant case, it is more convenient 
(see text) to work with block-sites and block-couplings shown 
here. 



Most of our discussion will use this pure imaginary 
form. In the absence of both SR-invariance and T- 
invariance, the different matrix elements of Him take on 
roughly independent imaginary values (apart from the 
requirements imposed by hermiticity). Thus, we have a 
general pure imaginary random hopping (ImRH) prob- 
lem on a doubled lattice. Each original fermion orbital 
a = is "represented" by two "copies" la and Ila; 

the energy of an orbital a is represented by an imaginary 
hopping amplitude between la and Ha, while hopping 
and pairing amplitudes between two orbitals a and f3 are 
represented by imaginary hopping amplitudes between 
the two pairs {la, Ila} and {1/3,11/3} (see Fig. ||a). 

In the presence of T-invariance, h and A satisfy 
<jyh*<ry = h and <yy/S*(Ty = A, where cTy = ay ® Ijv 
and CTy acts on the spin degree of freedom at each site. 
The corresponding restrictions on Him are best stated as 
follows: Group the different spin states on the same lat- 
tice site i and with the same copy index K into a block 
pair {Ki '\,Ki I}. T-invariance then implies that there 
are no internal couplings within such blocks. Moreover, 
the couplings between two blocks with the same copy 
index have a form 



f ia 


ib ^ 


— ir 


' coa0 


sin 
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, sin 6 


— cos 



(5) 



form 



with real c and d. Simply relabel- 



with real a and b {r — \/ a^ + 6^), while the couplings 
between two blocks with different copy indices have a 
i c —id 
id ic 

ing the spin states of copy II brings all block-couplings 
to the same form Eq. and we will find it conve- 
nient to discuss the T-invariant case using this picture 
of block-sites connected by block-couplings (see Fig. 
where each fat arrow represents the corresponding block 
coupling Eq. (^ with independent a and b, and differ- 
ent arrows are roughly independent of each other). We 
also use another representation of H in the T-invariant 
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case: This is obtainedEJ by performing a particle/hole 
and spin-up/down mixing unitary transformation Ur = 
^(l2(XiT2 + (Ty(8)rj:)(8)lAr (here, I2 and the PauU matrix 
ay act on the spin label, while the Pauli matrices r act 
on the particle- hole label): 

Tir thus has a bipartite complex hopping form, which 
will prove useful for our transfer matrix analysis (note 
however that strong correlations among matrix elements 
of Tir limit its usefulness in other contexts) . 

Finally, in the spinless case with no T-invariance, we 
simply drop the spin label altogether and do not have 
any constraints on the bonds of the corresponding ImRH 
problem. [Here and elsewhere in this article, we use 
the designation 'spinless' to also refer to all situations 
in which the spin label of the electron can be dropped — 
such as spin-polarized triplet superconductors. Also, in 
the rest of this article, we do not consider specifically the 
spinless case with T-invariance, but only mention the rek 
evant results as we go along. This case is very specialcl 
and maps onto a class of bipartite hoppinaM-oblems that 
have been studied both in onetil and twdla dimensions. 
In one dimension, our RG approach yields results con- 
sistent with what is known from Ref. ^ and can pro- 
vide more details about the low-energy properties of any 
particular system in this universality class, while some 
two-dimensional systems in this universality class will be 
discussed separately.tj] 

In this paper, we consider three one-dimensional sys- 
tems: the spinless superconductor without T-invariance, 
the 'spinful' T-invariant superconductor, and the 'spin- 
ful' superconductor without T-invariance. In all three 
examples, the bulk of our discussion is for the case with 
a single transverse mode (channel) present. However, 
as will become clear from our analysis of the 'spinful' 
superconductor without T-invariance, our basic conclu- 
sions regarding the low-energy physics apply equally well 
for the multi-channel case of all three problems. 

We model these systems by the appropriate BdG 
Hamiltonians with only nearest neighbor hopping and 
pairing amplitudes, in addition to same-lattice-site terms 
(the restriction to such nearest neighbor models is not at 
all crucial for any of our conclusions). The ImRH prob- 
lem corresponding to the single-channel spinless case is 
a two-leg ladder with all couplings pure-imaginary and 
roughly independent, but no other restrictions — in par- 
ticular, couplings along the rungs of the ladder are al- 
lowed. Such 'rung-couplings' (which we will sometimes 
refer to as 'vertical' couplings — see Fig. |l]a) correspond 
to any on-site terms in the original lattice BdG Hamil- 
tonian; in the spinless case, these can only be random 
potential terms, but more generally, one can also have 
"s-wave" pairing amplitudes and spin-flip scattering po- 
tentials (the latter due to frozen magnetic impurities). 



For the single-channel 'spinful' T-invariant case, the cor- 
responding ImRH problem is an analogous two-leg ladder 
of block-sites with pure imaginary block-couplings that 
are roughly independent of one another — in other words, 
we have a four-leg ladder with this special block struc- 
ture. Finally, the single-channel 'spinful' case without 
T-invariance is a pure-imaginary four-leg ladder with no 
other restrictions. These ladder problems are related to, 
the bipartite random hopping (RH) ladder problemsllfl 
in which rung couplings are disallowed but the other 
hopping amplitudes do not have to be pure imaginary — 
as mentioned earlier, these are in the same universality 
class as spinless superconducting wires with T-invariance. 
[The connection between this special class of models with 
sublattice symmetry and the more general systems that 
we study here is in fact surprisingly close, as will become 
apparent from our transfer matrix and RG analyses.] 

Finally, it is worth emphasizing at this point that our 
focus in all these cases is on the low-energy physics: The 
systems mentioned above will all differ from each other 
at higher energy scales, particularly in the initial 'dif- 
fusive' crossover regime that would be present in any 
weakly-disordered problem. However, at energies below 
this crossover scale, the effective value of disorder be- 
comes large and all these systems can be described by a 
unified physical picture. 

IV. STRONG-RANDOMNESS RG APPROACH 

In order to go beyond the heuristic ideas of Sec. we 
need a controlled approach that works in situations with 
strong Griffiths effects (such that z, defined by the low- 
energy behavior p(e) ~ g^i+i/^ of the density of states, 
is large). Such situations are expected to correspond to 
large values of effective disorder in the low-energy limit 
(in fact, a simple scaling argument indicates that the 
width of the distribution of the logarithms of the effective 
couplings is expected to be of order z in the low-energy 
limit). We are thus led to formulate a strong-disorder 
RG approach to this problem. 

Consider the ImRH Hamiltonian Him = X]y ^jjN)OI 
with tji = t*j — —tij. The eigenstates of ?iim occur in 
pairs with energies ±e, and the strong-randomness RG 
proceeds by eliminating, at each step, such a pair of states 
with energies at the top and at the bottom of the band: 
One finds the largest (in absolute value) coupling in the 
system, say ti2 connecting sites 1 and 2; this defines the 
bandwidth of the problem 2S7 = 2max{|ty|}. If the dis- 
tribution of the couplings is broad, the eigenfunctions of 
the two-site problem 7iiin[l,2] will be good approxima- 
tions to the eigenstates with energies ifi, since the cou- 
plings tij and t2j of the pair to the rest of the system will 
typically be much smaller. These couplings can then be 
treated perturbatively, and eliminating the high-energy 
states living on the sites 1 and 2 gives us the following 
effective couplings between the remaining sites: 
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) ^hj - ti2{t2l) ^tlj- 



(7) 



Clearly, the renormalized Hamiltonian Tiim again corre- 
sponds to a pure imaginary hopping problem, but with 
two fewer sit-in particular, the matrix Him has no di- 
agonal ternis,c3 in = 0. 

Some remarks on the proposed RG approach are in 
order here. Note that our RG rule Eq. (0) is an exact 
transformation for the zero-energy wavefunction, and as 
such provides information on the zero-energy localization 
properties. From a numerical point of view, it can be 
viewed as a construction of the zero-energy wavefunction 
in an a priori stable manner. Moreover, this transforma- 
tion can also be viewed as an approximate but accurate 
scheme for evaluating the 'Sturm sequence', i.e., the inte- 
grated density at very low energies, with the approxima- 
tions involved being well-controlled when the low-energy 
effective couplings are broadly distributed. Lastly, note 
that in the limit of strong disorder and in the absence of 
strong correlations between different couplings, the right 
hand-side of Eq. (|^) can be replaced with the "max" (in 
absolute value) of the three terms. 

So far, we have ignored the restrictions that would 
be imposed on the couplings of the ImRH problem in 
a 'spinful' (i.e, when both spin species need to be con- 
sidered) T-invariant situation. In this case, it is more 
natural to work with the block-sites and block-couplings 
defined earlier (Fig. |l]b). To begin with, note that the 
eigenstates now occur in doubly degenerate pairs with 
energies ±e. Thinking in terms of blocks automatically 
incorporates this degeneracy, since the eigenstates of a 
two-block-site problem Ht[1, 2] come in such doubly de- 
generate pairs with energies ±ri2 — ±\/ai2 + ^12- 
RG approach now eliminates four states at each step, 
two at the top and two at the bottom of the band with 
energies ±£7, where = maxlry}. The resulting effec- 
tive block-couplings among the remaining block-sites are 
again given by Eq. (|^) , but now each is a 2 x 2 ma- 
trix of the form (0). The effective problem is again an 
ImRH problem in the same block-form, and no block- 
diagonal terms are generated. For the RG rule Eq. (|^), 
the flows of bond-energies r and bond-angles 9 do not sep- 
arate. However, for strong disorder and in the absence 
of strong correlations (i.e., roughly, when the "max" RG 
rule applies) the energy variables r flow exactly as in the 
ImRH representation of a spinless superconductor with- 
out T-invariance, while the angle variables simply ran- 
domize. Thus, we expect essentially identical results for 
this 'spinful' problem with T-invariance and the corre- 
sponding spinless problem without T-invariance when- 
ever both flow to strong disorder sufficiently rapidly. 



V. CRITICAL POINTS AND GRIFFITHS 
EFFECTS IN ID: TRANSFER MATRIX 
ANALYSIS 

Before we go on to our more detailed RG analysis of 
the low-energy properties, it is useful to have a picture 
of the phase diagram in each case. We use an elemen- 
tary transfer matrix analysis to develop such a picture 
in terms of the zero-energy localization properties of the 
system. The goal here is to show by direct means that 
generically all our models are localized, and to demon- 
strate that there are critical points representing transi- 
tions between distinct localized phases; of course, these 
critical points can only be accessed by fine-tuning the 
disorder distributions. 

In general terms, we are looking at the zero-energy 
localization properties of M-leg ladder systems governed 
by a Schroedinger equation 



(8) 



where in and Un are M x M matrices, and the wave- 
function tpn is an M-dimensional vector defined at each 
rung n. We find it convenient to work with the following 
transfer matrix formulation: 



a M.tt.^Vi 



this defines elementary transfer matrix r„. 



1. The bi-partite ladder problems: 

To begin with, consider the bipartite problem men- 
tioned above, in which one has M coupled chains with 
no rung couplings, = 0. This system was studied by 
Brouwcr et. a/.jlfl who found that for both real or com- 
plex hopping amplitudes there are M localized phases 
separated by M dimerization driven delocalized critical 
points; each critical point exhibits a strong Dyson singu- 
larity in the density of states, p{e) l/£\ ln"^e|. Here, we 
rederive by completely elementary means the existence 
of M delocalized critical points, and also characterize the 
M + 1 distinct localized phases; the ideas introduced in 
the process will generalize naturally to our superconduc- 
tor systems. 

Being bipartite, the two sublattices decouple at e = 0: 



T,i+iT„ — 



_f-i ft 








(10) 



We are thus led to study the Lyapunov spectrum of the 
matrix products ^^^.^t 



2k ^2k-l 



t2^t[ and p2f.t 



2k''2k-l 



. 121-1 ■ 



At zero dimerization (t2fe and ^2/0+1 distributed identi- 
cally) the Lyapunov spectra of both products (consisting 
of M distinct Lyapunov exponents in the general case) 
are identical and symmetric around zero, and the full 
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Lyapunov spectrum (of the full transfer matrix product) 
is doubly degenerate. Thus, for even M the smallest (in 
absolute value) Lyapunov exponent is non-zero and the 
zero-energy modes are localized, while for odd M there 
is always a zero Lyapunov exponent. [The actual values 
of these exponents are not of much concern here, only 
the fact that they are all distinct.] Now, consider adding 
dimerization by simply multiplying every odd bond i2fe+i 
by a scale factor e ^ and every even bond t2fc by a fac- 
tor . Clearly, the whole Lyapunov spectrum for one 
reduced (sublattice) problem is shifted rigidly by 5 (and 
by the exactly opposite amount for the other sublattice). 
Thus, as we scan 5 from — oo to -l-oo, there will be ex- 
actly M points where two of these exponents of the full 
transfer matrix cross zero. For the given sublattice, if 
we label each non-critical region by {k,M — k) with k 
growing modes and M — k decaying modes, these crit- 
ical points represent consecutive "delocalization" tran- 
sitions (k, M — k) ^ {k + 1, M — k — 1). Of course, this 
can also be restated in terms of the number of zero-energy 
states localized at each end in a finite odd-length chain 
with free boundary conditions, and provides a 'topolog- 
ical' distinction between the different localized phases; 
each transition corresponds to a single zero-energy state 
becoming delocalized and migrating from one end of the 
chain to the other. 

In a bulk system, at any of these critical points, one 
has two zero-energy delocalized Lyapunov modes (linear 
combinations of which can roughly be interpreted as a 
'left-moving' and a 'right-moving' slow modes). Since 
any description of the low-energy properties of the crit- 
ical system in terms of such 'slow' modes has to respect 
the bipartite nature of the original problem, it is natural 
to expect that such a low-energy effective theory is in 
the same uniitersality class as the bipartite single chain 
RH problem.Ej This provides a clear rationale for the 
critical low-energy density of states to be of the Dyson 
form. Similarly, we expect that the localized phases in 
the vicinity of any critical point look, at low energies, 
like the dimerized Griffiths phases obtained by introduc- 
ing a small amount of dimerization in the single chain 
problem — in particular, we expect Griffiths singularities 
in the density of states consistent with the results of 
Ref. 0. 

2. Single- channel spmless superconductor without 
T-tnvariance: 

Returning to the dirty superconductor problems, con- 
sider first the spinlcss fermion system with no T- 
invariance in the ImRH language. In this case, the rung 
coupling is of the form Un = fJ-n^y with some real fi„ (ay 
acts on the copy label K), while the hopping term is a 
general pure imaginary 2x2 matrix in . Because of the 
identity aya~^ay — oF j deta valid for any 2x2 matrix 
a, the zero-energy transfer matrices "decouple" : 



lo oy) U„f-V, ; 1,0 aj' ^^^^ 

where r„ = dett„. Thus, Lyapunov exponents of the 
product of T„ are given by a "superposition" of the Lya- 
punov exponents of the product of \/\Tn\in^cry and of 

the product of ( ^"/\^^^ ^he former 

Vsign(T„)Vr„| / 

product is very similar to the product ^ki2ki2k-i stud- 
ied earlier (note, however, the particular "normalization" 
used here) ; the Lyapunov spectrum consists of two expo- 
nents ±7t, with 7t of order one. The latter product is 
essentially the transfer matrix product at e = for a Id 
Anderson problem with random energies ±/^„ (depend- 
ing on the sign of t„) and random hopping amplitudes 
^ I r„ I ; the two corresponding exponents are ±7^. The 
full Lyapunov spectrum thus consists of the four expo- 
nents ±7t ± 7^ . As we increase the strength of the rung 
couplings /i„ (e.g., by increasing the root mean square 

strength R = i?^ = with JI = kept fixed) from 

to cx), 7^ also increases from to 00. Thus, at some 
critical strength R = Rc oi the rung couplings, 7^ will 
equal 7t, and two Lyapunov exponents will be zero cor- 
responding to an isolated delocalized critical point along 
the R axis. 

Now, in the superconductor problem, there are no 
eigenstates at precisely zero energy in any system — this is 
due to the presence of the rung couplings. Nevertheless, 
there may be states with exponentially small (in system 
size) energy, and the localized phases on either side of 
the critical point can again be characterized in terms of 
such almost-zero-energy end-states. Consider again an 
open odd- length chain. When yUn = (and i5 = 0), there 
are two zero-energy states, one at each end of the chain. 
Turning on the rung couplings enables the two to mix, 
but since they are separated by the entire length of the 
chain, the splitting is exponentially small in the length of 
the chain — as long as i? < Rc- Thus, in this phase, there 
will be two such (essentially) zero-energy end-states. In 
terms of the quasiparticle spectrum of the original su- 
perconductor Hsc, there is a single quasiparticle state 
with an exponentially small energy and a wavefunction 
with weight only at the two ends of the chain. On the 
other hand, in the phase with R > R^ there are no such 
end-states with nearly zero energy, as may be argued 
by starting with the limit of large /x„. Thus, we have 
two different localized phases distinguished by this topo- 
logical property; if we think in terms of a more general 
'dimerization'-'rung coupling' (S — R) parameter plane, 
we have a phase diagram shown schematically in Fig |^. 
Moreover, we again expect the low-energy properties in 
the vicinity of the transition between the phases to be 
in the universality class of the single chain bipartite RH 
problem in the vicinity of its zero-dimerization critical 
point. In particular, we again expect Griffiths singular- 
ities in the density of states of either localized phase in 
the vicinity of the critical point — our expectations will 
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FIG. 2. Schematic phase diagram in the S ~ R plane for 
both the general spinless case without T-invariance, and the 
'spinful' model with T-invariance. We use R to denote some 
measure of the strength of the rung couplings which corre- 
spond to the random on-site terms of the original supercon- 
ductor problem (e.g., for a symmetric distribution of random 
rung couplings, we can define _R to be a root mean square 
of this distribution). The vicinity of the phase boundary is 
expected to exhibit strong Griffiths effects. 

be borne out by the more detailed RG analysis in the 
next section. 



3. Single- channel spinful T-invariant superconductor: 

Consider now the 'spinful' system with T-invariance 
(but no SR-invariance). In the bipartite representation, 

Eq. ^ , the rung coupling is u„ = ( f Wncry \ .^-^j^ 
some c- number Wn, while the hopping term has a form 
£„ = ^ '^Q^ with hn = —a* but otherwise gen- 
eral complex 2x2 matrix a„. From the transfer matrix, 
Eq. (^), the sublattice decoupling at zero energy is seen 
immediately; on one sublattice, we need only consider 
the product of 

Lyapunov exponents of this product together with 
their negatives (from the other sublattice) form the 
full Lyapunov spectrum of the system. But this 
product further "decouples" precisely as in Eq. ( [ll| ) 
into the product of ^J\an\a~^ay and the product of 

fwj^ l/^Vherea„ = deta„. Lyapunov 

exponents of the former product are ±7^ with some 70 of 
order one. However, the Id Anderson localization prob- 
lem corresponding to the latter product is non-hcrmitian: 
while the hopping amplitudes can be chosen real and 
equal to A/|a„|, the on-site energies |wn|e*'^" are com- 
plex, with the phases ipn having contributions from the 
phases of both w and a. Nevertheless, the Lyapunov 
spectrum still consists of two exponents ±7tu. The spec- 
trum of the product of the Q„ is thus ±7^ ± 7^,, and 



the other sublattice merely duplicates this to make the 
full Lyapunov spectrum doubly degenerate. Now, all we 
need to know is that for small values of rung-couplings 
the corresponding 7^ is also small, while for large values 
of rung-couplings 7^ is large. Then, as in the spinless 
case, there is an isolated delocalized critical point along 

the R = Ryj = \/'vip' axis for some critical strength Rc of 
the rung coupling terms (at which —la)- Note that at 
this critical point, a total of four Lyapunov exponents will 
simultaneously cross zero; the corresponding two pairs of 
critical modes are related to each other by T-invariance. 
Again, the phase with R < Rc in characterized by the 
presence of end-states with exponentially small energies. 
The only difference from the spinless case is that there are 
now four of them — this corresponds, in terms of Hsc, to 
two T-symmetry related quasiparticle states, each with 
an exponentially small energy and a wavefunction having 
weight only at the two ends of the chain. The phase with 
R > Rc again has no such nearly-zero-energy end-states. 
In the full S-R plane, we thus have the schematic phase 
diagram shown in Fig ^j. Again, the critical point and 
the phases in its vicinity are expected to look, at low en- 
ergies, like those in the single chain RH problem, with 
an additional degeneracy introduced by T-invariance. A 
direct numerical check of the Dyson form for the criti- 
cal density of states in this 'spinful' T-invariant case is 
shown in Fig. ||, while the RG results of the next section 
confirm the underlying physical picture in detail. 

Parenthetically, we also note that these Lyapunov ex- 
ponent crossing arguments imply that in all cases consid- 
ered so far, the inverse of the smallest Lyapunov exponent 
diverges as \g — gc\~^ where g is some tuning parameter 
that drives the system through the critical point at gc- 
This implies that the exponent v for the typical localiza- 
tion length is J/ = 1 at all these transitions. 

4- Single- channel spinful superconductor without 
T-invariance: 

Finally, we turn to the 'spinful' case with neither T- 
invariance nor SR-invariance (again, with only a sin- 
gle transverse channel). [Note that a quasi-lc? spin- 
less fermion superconductor with two transverse channels 
would also be modeled by the same BdG problem, with 
a somewhat different interpretation for the various cou- 
plings. Thus, analysis of this case is of value in demon- 
strating that our basic conclusions are not special to the 
single-channel case in any of the problems we consider.] 

In this case, we have been unable to come up with any 
simple decoupling scheme that maps the corresponding 
transfer matrix to that of some problem with sublattice 
symmetry — the exact mapping we have used earlier is 
thus special to the two single-channel cases considered 
above. However, when the rung couplings are all zero, we 
do know that the corresponding bipartite four-leg ladder 
has five phases as we scan the dimerization parameter 
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FIG. 3. Numerical check of the Dyson singularity in the 
'spinful' T-invariant case. Number of states with energies be- 
tween and E for an open chain of length L = 10^, averaged 
over 5 disorder realizations, is plotted vs the log-energy scale 
r_B = ln(l/i5). Independent inter-site couplings (real and 
imaginary parts, of independent hopping and pairing ampli- 
tudes) are chosen from a uniform distribution over [—1,1], 
while independent on-site couplings (real chemical potential 
and real singlet pairing amplitude) are chosen from a uni- 
form distribution over [-W, W]. Critical Wc = 1.630 ± 0.001 
was found accurately from numerical transfer matrix analy- 
sis. At this point, we clearly have A^r ~ l/F^ (which is shown 
displaced from the data for clarity). We also show several 
off-critical points and give rough estimates of the correspond- 
ing dynamical exponents z from the Griffiths fits A''^; ~ iJ'*/^ 
over 4 < Ai5 < 1000. 

5 — these are labeled (fc,4 — k) with fc = 0, 1,2,3,4 cor- 
responding to 4 — fc zero-energy states localized at one 
end and k states localized at the other end for a finite 
ladder with an odd number of rungs and free boundary 
conditions. In particular, in the vicinity of (5 = one has 
two such states at each end. Now, turning on some weak 
rung couplings allows these two states (at each end) to 
mix amongst themselves, and there will thus be no states 
with exponentially small energies in this regime. On the 
other hand, the same is clearly true for very large values 
of the rung couplings. Thus, the phases obtained in ei- 
ther case in the vicinity oi S — are expected to 'look' 
the same (this is made more precise later using our RG 
approach) . 

The question then arises: Is there a single intervening 
phase (at (5 = 0) which is topologically different? This 
would result in two phase transitions as we scan the mag- 
nitude of the rung couplings. Another simple possibility 
is that there is no transition at all as a function of in- 
creasing rung coupling. Of course, there can also be other 
more complicated scenarios (some possibilities for the full 
phase diagram are sketched in Fig. Moreover, since 
there are several independent rung couplings correspond- 
ing to each physical lattice site, there are many different 
ways of 'increasing the rung coupling' and the possible 



FIG. 4. Possible phase diagrams in the 5 — R plane for the 
'spinful' case without T-symmetry. We use R to represent 
some measure of the strength of the rung-couplings. On the 
diagrams, we denote the phases with zero-energy end-states 
by 1 and the phases with no such states by 0. A and B rep- 
resent two simple possibilities, while a more complicated case 
with a multi-critical point is shown in C. Our numerical trans- 
fer matrix studies tentatively suggest that in the most 'ran- 
dom' such superconductor (i.e., with all couplings present and 
independent) the phase diagram is very nearly that of the 
panel C. In many other cases (e.g., when we allow only on-site 
chemical potentials in the original superconductor problem, 
but completely general inter-site couplings), we observe the 
phase diagram A. We have not found realizations that would 
clearly exhibit the phase diagram B; however, we do not have 
arguments that would rule out this possibility. 

phases and transitions encountered along the way will 
most likely depend on how we scan. 

To get a more detailed picture, we performed an ex- 
tensive numerical transfer matrix analysis (to obtain ac- 
curate results, we used the technique of Ref . ^) , as well 
as exact diagonalization studies, for particular choices 
of scan. In one choice of scan, we include all possible 
inter-rung couplings, but allow only those (intra-)rung 
couplings that correspond to random on-site chemical 
potentials in the original superconductor problem. In 
this case, we can clearly delineate the phase boundaries 
to conclude that we have a phase diagram of the type 
shown in Fig. with no transition along the R axis 
at (5 = (with R now a measure of the strength of the 
on-site potentials). If we scan across the phase bound- 
ary at a fixed non-zero value of dimerization (so that we 
start in the phase with zero-energy end-states and leave 
it by increasing the mean-square strength of the random 
potentials), we find that the typical localization length 
defined by the inverse Lyapunov exponent again diverges 
with an exponent v — 1, and the critical density of states 
is again of the Dyson form. 

In another scan, we take the most random such system 
(i.e., with all allowed couplings present and independent) 
and boost the strength of all the (independent) rung cou- 
plings by the same amount on average, while keeping 
distribution of the inter-rung couplings fixed. As we in- 
crease the typical magnitudes of the rung couplings from 



zero in this manner, we again see strong Griffiths singu- 
larities in the density of states developing, and, possibly, 
critical behavior. However, near such a tentative criti- 
cal point the two Lyapunov exponents that come close 
to zero seem to have a mutual "level repulsion" and the 
Lyapunov spectrum seems to exhibit the analog of an 
"avoided level crossing" ; this results in huge localization 
lengths and strong near-critical behavior in large regions 
around the tentative critical point. A more detailed in- 
vestigation away from i5 = yields a phase diagram of 
the type shown in Fig. Thus, to within our numerical 
accuracy, there seems to be a multicritical point in the 
5 — R phase diagram at S — 0. However, we can not 
exclude the possibility that we are seeing a case with no 
transition along the i?-axis (but phase-boundaries com- 
ing very close to this axis) or an almost-degenerate case 
with two very closely spaced transitions. Crossing the 
phase boundaries away from the putative multi-critical 
point again gives a localization length exponent of = 1 . 
However, we are unable to make any reliable statements 
in the vicinity of the multicritical point. 

We have scanned along several other directions in the 
parameter space (corresponding to different interpreta- 
tions of the rung-coupling parameter i?), but have not 
clearly seen two distinct transitions as in Fig. ^B; how- 
ever, behavior of the type shown in Fig. ^ is the most 
commoiijHowever, note that the weak-disorder analysis 
of Rcf. ^ did find critical behavior in the conductance 
at S — with inter-rung and intra-rung couplings iden- 
tically distributed and chosen from a Gaussian distribu- 
tion. Moreover, their result is consistent with our predic- 
tions for critical behaviour at 'ordinary' critical points, 
as opposed to multicritical points.] 

These numerical results thus confirm our suspicion 
that the structure of the full phase diagram in this gen- 
eral 'spinful' case (or in multi-channel versions of all the 
cases) is quite complicated. The specific phase diagram 
obtained by tuning the parameters of a particular physi- 
cal system (in which some subset of the allowed couplings 
may be identically zero) can thus be very different from 
case to case (a similar observation in related 2d mod- 
els has been made in Refs. ^ and |2^). However, it is 
clear that there can be, in general, two kinds of localized 
phases. Moreover, whenever both phases are present in 
the phase diagram of a particular system, we again ex- 
pect (analogous to the single-channel cases considered 
above) the system in the vicinity of the phase bound- 
ary to 'look' at low-energies like a single-chain random 
hopping problem with weak dimerization (note that the 
numerical estimate, v ^ 1, that we obtain away from 
any special 'multicritical' points supports this picture, 
and our RG results provide further confirmation of the 
same) . 

This is as far as we can go with an elementary transfer 
matrix analysis. For a more detailed characterization 
of the low-energy properties, we now turn to the strong 
disorder RG analysis. 



VI. CRITICAL POINTS AND GRIFFITHS 
EFFECTS IN ID: RG ANALYSIS 

To test the above picture in detail, we have imple- 
mented the RG numerically in all three single-channel 
cases. Since the single-channel spinful problem without 
T-invariance is essentially identical to the correspond- 
ing two-channel spinless problem, this last example also 
serves to establish that our conclusions do not depend on 
any special properties of the single-channel case. In the 
single-channel spinless case without T-symmetry (and 
the corresponding spinful case with T-invariance) , we fo- 
cus mainly on the vicinity of the phase transition at 6 = 
(see Fig. |]). In the 'spinful' case without T-invariance, 
we consider the two different realizations described in 
the previous section — these have phase diagrams of the 
types shown in Fig. |^A (in which we scan across the phase 
boundary at non-zero dimerization so that we leave the 
phase having end-states by increasing the strength of the 
onsite potentials) and in Fig. ^ (in this case, we focus 
on the immediate vicinity of the putative multicritical 
point). [For completeness, we have also studied the tran- 
sitions as a function of di merization in the bipartite lad- 
der problems of Sec. |V l[ ] 

For the S — spinless case, (and the 'spinful' T- 
invariant case) we consider systems with lengths as large 
as L 2 X 10^. The initial conditions used have ran- 
dom (inter-rung) hoppings chosen from a uniform distri- 
bution over [—1,1], and random (intra-)rung couplings 
u chosen from a symmetric (with either sign equally 
probable) power law distribution P{u) — 
|u| < 1. Note that the RG transformation (|^) can be for- 
mulated entirely in terms of the imaginary parts of the 
couplings, and this is the language that we use here. For 
the 'spinful' case without T-symmetry, we are restricted 
to L < 5 X 10''^. When we scan at finite dimerization, 
we enforce this dimerization by choosing the even and 
odd inter-rung couplings from the power law distribu- 
tion P{u), but with different fixed values gcvcn 7^ 5odd 
for the power law exponents. The strength of the rung 
couplings is again tuned by varying the corresponding 
power-law exponent (?rung- We use the "full" RG rule (|^) 
rather than its less accurate "max" version since we are 
primarily interested in testing our physical picture for the 
low-energy properties starting with a system with mod- 
erate values of the bare disorder. 

Apart from the immediate vicinity of the putative mul- 
ticritical point of Fig. ^ (which we comment upon sepa- 
rately), our results are equally reliable and essentially 
identical in all the cases studied. In the interests of 
brevity, below we display in detail only the results ob- 
tained in the spinless case at (5 = (see Fig. ||). 

We search for the critical point by looking at the frac- 
tion of "isolated" sites among the remaining sites in 
the system: If there are N sites left, we find all the 
sites that are "covered" by the N/2 strongest remaining 
bonds. Roughly speaking, the fate of these covered sites 
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FIG. 5. Fraction of isolated sites among the remaining sites 
for the two-leg ImRH ladder representing spinless fermion su- 
perconducting wire. Critical gc ~ 0.95 is clearly identified. 
Inset shows the fraction of cluster decimations (see text for 
details). 

is clear — they will be frozen dynamically at rougly this 
log-energy scale. The rest of the sites are still dynami- 
cally free at this scale, and we call them "isolated". In 
a localized phase, this fraction quickly approaches zero. 
On the other hand, at a critical point, we expect that new 
couplings which contribute, upon their subsequent deci- 
mation at a lower energy, to the density of states at that 
lower scale, are formed continually over all energy scales 
in a scale-invariant way. The fraction of isolated sites 
at criticality is therefore expected to saturate to some 
constant at low energies. In the case of a single critical 
RH chain, a quarter of the remaining sites is notionally 
isolated at each stage of the RG: each site has a bond 
to the right and to the left, and the bond strengths are 
uncorrelated; since each bond has a 50% chance of being 
"weak", the site is "isolated" with probability 1/4. This 
gives us a very direct test for the location and nature of 
the critical point. 

Results of such a search for our spinless case are shown 
on Fig. H We clearly identify a critical gc ~ 0.95, and 
note that the fraction of isolated sites at this critical 
point remains essentially 0.25, which is evidence that the 
low-energy behavior is that of a single one-dimensional 
'backbone' that goes critical. Note that it is also pos- 
sible to further probe the geometry associated with the 
low-energy theory, as is done, e.g., below in our "order- 
parameter" studies, and even more exhaustively by look- 
ing in detail at probability distributions of the various 
couplings (which we have not pursued fully — the RG re- 
sults shown here already confirm our basic picture of the 
low-energy physics) . 

At the critical point, the number of sites remaining at 
the log-energy scale F = ln{flo/fl) is iVr ^ 1/F^, as can 
be readily seen from Fig. ^j. Since iVr is essentially the 
integrated density of states, this is consistent with the 



FIG. 6. Number of remaining sites A'^ vs log-energy scale 
r. Also shown is the fit of the critical Nr to the form A/V'^. 

Dyson p(e) ^ l/e|ln'^e|. Moreover, at g — gc, the dis- 
tributions of all couplings become broader and broader 
on the logarithmic scale, with the characteristic widths 
growing linearly with F. This is shown on Fig. ^. In con- 
junction with our result for the fraction of isolated sites 
at low energies at criticality, this indicates that the crit- 
ical point is in the same universality class as the single 
chain RH problem. [Note that since the system effec- 
tively reduces to a single RH chain 'back-bone' in the 
low-energy limit, we expecIO the critical average ther- 
mal conductance Kc{L) of the original superconductor to 
scale as kBT/^/h, where L is the length of the wire — 
this is consistent with the weak disorder result of Ref. ||. 
As mentioned earlier, this leads us to believe that their 
analysis was performed only at such a critical point and 
does not represent the generic behavior of the system.] 

Consider now the two different localized Griffiths 
phases. As expected, we find A^r ^ e^^^^ with a 
continuously varying dynamical exponent z which di- 
verges as z ~ \g — gc\~^ (to within our numerical ac- 
curacy) as we tune across the critical point. Within our 
strong-disorder approach, the two Griffiths phases are 
distinguished by the character of the corresponding RG- 
generated dimer patterns. To discuss this intuitive dis- 
tinction more precisely, we compare such RG-generated 
dimer covers against some fixed reference dimer cover. A 
natural choice of the reference cover for the spinless prob- 
lem at hand is a "vertical" cover with reference dimers 
covering the rungs of the ladder, i.e., joining the two 
copies of each of the original fcrmions of Hsc- We use 
this specific reference cover in our numerical studies pre- 
sented below; the following discussion, however, is fairly 
general. [We use the same 'vertical' reference cover in 
the other two superconductor problems we study with 
the RG.] 

We call two sites connected by such a reference bond a 
cluster, and the corresponding coupling between the two 
sites a field on the cluster. Couplings that connect sites in 
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FIG. 7. At the critical point (determined from Fig. |5|) dif- 
ferent measures of the widths of the tog-couphng distributions 
aU scale linearly with F. We plotted the A'^/2-width of the dis- 
tribution of all bonds, the 'decimation error' width (defined 
as the average logarithm of the ratio of the decimated bond 
to the strongest nearby bond), and also the A'^/4- widths of the 
interaction and field distributions with respect to the vertical 
reference cover (see text for details). 



different clusters are called interactions between the clus- 
ters. This terminology is borrowed from the Id random 
transverse field Ising model (RTFIM), but we empha- 
size that the correspondence is not exact although we do 
expect that the critical behavior, characterized with re- 
spect to a well-chosen cover, is essentially that of the RT- 
FIM (we expect this to be true because a similar analysis 
for the dimerized single chain random hopping problem 
using the natural reference cover consisting of alternate 
bonds gives an exact mapping to the Id RTFIM). When 
a coupling connecting two sites in the same cluster is 
decimated (field decimation), the corresponding cluster 
is killed. When a coupling connecting two sites in dif- 
ferent clusters is decimated (interaction decimation) , the 
two clusters are joined into one new cluster, which is now 
specified by the two other (remaining) end-points of the 
two original clusters. The number of original sites that 
belong to a cluster defines its 'magnetic moment' and it 
then makes sense to talk of a magnetization density m for 
the system. Pictorially, a reference dimer cover specifies 
some connectivity rules for the RG-formed dimers, with 
natural notions of clusters and percolation with respect 
to such connectivity rules, while the RG rules prescribe 
the dynamics (as a function of increasing F) obeyed by 
the clusters. Such connectivity properties can be used to 
distinguish between different RG- generated dinter covers, 
and, thus, to characterize the different phases. E3 Within 
the strong-disorder RG, this is the distinction that cap- 
tures the different 'topological' character of the two Grif- 
fiths phases. 

Going back to our numerical RG studies, we find that, 
with respect to the vertical reference cover, in the phase 




FIG. 8. Average moment of the remaining clusters with 
respect to the vertical reference cover. Also shown is 
the fit at the critical point for the exponent 0. From 
this fit we obtain (j) « 1.59, which can be compared with 
the corresponding exponent for the quantum Ising model, 
0RTFIM = (1 + V5)/2 ^ 1.62. 



which obtains for strong rung couplings (small values of 
g) there are only small clusters and no connectedness 
across the whole system, while for weak rung couplings 
(large values of g) there is an infinite percolating cluster 
that forms in the limit of large F. This development of 
topological order for g > g^ is characterized by an expo- 
nent (3 defined by the scaling of the average magetization 
density: m{r — s- oo) ^ {g — gc)^ ■ At criticality, the av- 
erage moment of the surviving clusters scales as fi ^ V^, 
which defijies the exponent 0; in complete analogy to the 
RTFIM ,E3 the exponent (3 for the topological order pa- 
rameter ('magnetization') can be obtained from via the 
scaling relation [3 = 2 — (j). Figure || shows our numerical 
result for the F dependence of the average moment of 
the surviving clusters at the critical point. The numeri- 
cal value obtained for the exponent (j) is very close to that 
of the Id RTFIM. The value /3 « 0.41 we infer using the 
scaling relation is then very close to the corresponding 
exact result for the single dimerized chain (which has a (3 
exactly equatJ to-the magnetization exponent (3— ■\/5)/2 
of the Id RTFIMEI). We also note that, similarly to the 
quantum Ising model, the critical point is the point of 
balance between the cluster fields and cluster interac- 
tions; this is shown in the inset of Fig. |^, and provides 
alternative means for identifying the critical point. This 
completes our RG description of the spinlcss case. 

Essentially the same results are obtained for the corre- 
sponding spinful problem with T-symmetry — this is con- 

We 



sistent with our general argument in Section IV. 
also repeated this analysis for the specific realization of 
the spinful problem without T-symmetry corresponding 
to Fig. All the results obtained for this case in the 
vicinity of the transition at fixed finite dimerization (as a 
function of increasing R) are essentially identical to the 
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results shown above for the spinless case. [Entirely analo- 
gous results are also obtained in all cases for the dimeriza- 
tion driven transitions of the bipartite ladder problems.] 
Thus, our general picture for the low-energy physics ap- 
pears to be validated by the RG results so long as we are 
not in the vicinity of any special multicritical points. 

Finally, a brief comment on the RG results in the vicin- 
ity of the putative multicritical point. Fig. As we scan 
R at (5 = 0, strong Griffiths effects again show up clearly 
over a wide region in the vicinity of the putative multi- 
critical point. Moreover, the phases at large and small R 
both look 'paramagnetic' in terms of clusters defined with 
respect to the vertical reference cover (and also many 
other reference covers). This is consistent with our argu- 
ments in the previous section. However, our RG results 
are also inconclusive when it comes to pinning down the 
structure of the phase diagram near this apparent mul- 
ticritical point — again the analysis is plagued by near- 
critical behavior over a wide region. The corresponding 
long crossovers do not allow us to make any definitive 
statements regarding either the presence or the univer- 
sality class of such multicritical points. This remains the 
principal unresolved question in our entire analysis of the 
one dimensional examples. 



VII. HOW ARE THE LOW-ENERGY STATES 
GENERATED? 

We now consider precisely how the states in the singu- 
lar low-energy tail of the density of states are generated, 
and identify the corresponding Griffiths regions in the 
ImRH language. As is already implicit in the above dis- 
cussion, such a pair of low-energy states is formed if there 
are two 'isolated' sites separated by a region in which all 
the other sites are locked into short-ranged (dimer) pairs. 
Griffiths effects are a consequence of such isolated sites 
being generated sufficiently often. An example of a pair 
of such sites is shown in Fig. ^. The isolated sites are 
very weakly coupled to each other — the coupling is of or- 
der e ^ e~'^', where I is the length of the region. In a 
disordered system, there is always a probability of order 
p', with somep < 1, of finding such a region — this gives a 
power-law contribution ~ el'^ipl/c to the integrated den- 
sity of states from such regions. Thus, we always expect 
a variable power-law (Griffiths) density of states in such 
random hopping problems in which no on-site energies 
(i.e., diagonal terms) are allowed; we conclude that Grif- 
fiths phases are generic. The specific example shown in 
the figure is expected to be relevant to Griffiths effects 
in the phase in which the rung couplings dominate on 
average. The regions to the left and to the right of the 
pair are intended to be a caricature of the typical re- 
gions in this regime, while the region in the middle is 
comparatively rare. Decreasing the strength of the rung 
couplings would increase the probability of having such 
regions, corresponding to the observed increase in the 
dynamical exponent as one approaches the transition to 



FIG. 9. Griffiths generation of low-energy states. 



the other phase. The critical point is then characterized 
by a proliferation of such Griffiths regions on all energy 
scales. 

The picture that emerges is thus very similar to that in 
our toy- model of Sec. || in which the low-energy states 
are closely associated with 'domain walls' between the 
two different gapped phases of the pure system. 

It is useful, at this point, to recast some of this in terms 
of the original superconductor Hamiltonian Hsc- This 
will give us a somewhat different perspective on the origin 
of these low-energy states. For simplicity, we restrict the 
discussion below to the spinless case. We first examine 
the basic RG rule Eq. (|^). Consider a two-site problem 

H[a, P] = eaclca + ^pc^pCfj + {tcl^c/i + Aclcj^ + h.c). 

This can be diagonalizeji by an appropriate Bogoliubov 
transformation to givecj 

H[a, P] = £+7^7+ + e_7i7_ + const. 



e± = 



+ epy + 4|A|2 ± V(ea - 6/3)2 + 4|i|^ 



From this solution, it is clear that states at both a and f3 
simultaneously contribute to either eigenstate only if ei- 
ther \t\ > \ea~efi\ or |A| > |eQ-fe/3|. When this happens, 
we can no longer think of the two sites in isolation and 
need to account for the resonance between them. Now, 
if |tp — |Ap ~ ea^/i, £- for such a resonance can be very 
small. In order to obtain a good low-energy description 
in such a situation, we would eliminate the high energy 
'-I-' state and keep only the ' — ' state by introducing a 
single effective site with site energy e_. In our ImRH 
RG language, this corresponds to a situation in which 
a single bond connecting sites on two different rungs is 
decimated, leaving behind one site on each rung; these 
two remaining sites are coupled by a weak bond precisely 
equal in magnitude to e_ . On the other hand, if there 
is no mixing of the states at a and /?, i.e., when one of 
the local potentials, say ea, dominates, we would elimi- 
nate the state a completely; in the ImRH RG this cor- 
responds to decimating the corresponding rung. Thus, 
our RG procedure either eliminates a full fermion state 
because it is frozen out by a strong local potential, or 
eliminates 'half of a fermion' from each of the two sites 
in resonance and recombines the remaining halves to cre- 
ate a new effective fermion with site energy equal to the 
new coupling introduced (our RG is thus really defined 
on the corresponding Majorana fermion states). 

Now, consider for concreteness the Griffiths phase in 
which the on-site potentials dominate. At low enough 
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energies, a description in terms of isolated effective sites 
(with negligible mixing between them) with some renor- 
malized distribution of effective site energies is clearly ap- 
propriate. However, to arrive at such a description, one 
has to first account for all the resonances at higher energy 
scales that arise from any anomalous regions in which 
hopping and pairing amplitudes are large compared to 
local potentials. As a specific example of such anoma- 
lies, consider the central region of Fig. |[ In the origi- 
nal superconductor language, this region corresponds to 
both A and t being relatively large and comparable in 
magnitude, in addition to having a somewhat definite re- 
lationship between their phases throughout this region. 
Eliminating all the resonances between the sites in this 
region finally gives a low-energy description in terms of 
a single effective site with an exponentially small energy. 
This effective site in the original superconductor language 
clearly corresponds, in the ImRH language, to the pair of 
isolated sites shown in Fig. |[ The ImRH RG thus pro- 
vides the natural language for capturing the important 
low-energy physics. Crudely speaking, the 'effective site- 
energies' in the original superconductor language corre- 
spond to the 'cluster field couplings' with respect to the 
natural vertical reference dimer cover in the ImRH RG 
language. Note also that an important ingredient of the 
physical picture that emerges is the spatial character of 
the quasiparticle states (in particular, note that wave- 
function of a low-energy quasiparticle in the phase with 
zero-energy end-states is split into two spatially sepa- 
rated pieces) . Such spatial information is also kept most 
naturally in the ImRH RG; in particular, the develop- 
ment of the 'topological order' is seen most naturally in 
this language. 

Finally, we can now go back and ask what is the pre- 
cise role played by the various symmetry restrictions. As 
discussed above, Griffiths effects in the ImRH RG lan- 
guage are associated with situations in which we repeat- 
edly eliminate only a subset of the couplings connect- 
ing two rungs (for instance, only one of the couplings in 
the spinless example above), thus splitting the original 
fermion states. Now, if some symmetry restrictions re- 
quire that some of the original couplings between rungs 
equal one another, then we would be forced to decimate 
several of them simultaneously — for example, we could be 
forced to eliminate only complete fermion states. Within 
the RG approach, this could then rule out the possibility 
of Griffiths effects. It is easy to see that the restrictions 
imposed by T-invariance are not enough for this to hap- 
pen. On the other hand, in the SR-invariant case, a 
simple analysis of the symmetry constraints for the cor- 
responding IrnRH problem shows that this is precisely 
what happenscj In the strong disorder limit, this, then, 
is the true significance of the absence or presence of SR- 
invariance for quasiparticles of a superconductor. 



VIII. DISCUSSION 

In this article, we have established the existence of 
strong Griffiths effects in one-dimensional superconduct- 
ing wires in which the quasiparticle spin is not a good 
quantum number. We associated these singularities with 
quasiparticle states that live on 'domain walls' between 
adjacent large regions in two different phases, one phase 
that supports zero modes localized near the ends of a 
finite system, and another that does not support such 
modes. 

An obvious question now arises: Do such effects exist 
in two or more dimensions? Thinking in terms of the 
ImRH RG, it does seem that such effects could exist in 
cases without SR-invariance, particularly in the vicinity 
of the thermal-metal to thermal-insulator transition, or 
the transition between the two topologically distinct in- 
sulating phases in cases without T-symmetry. However, 
it is of course not a priori clear if 'isolated sites' would be 
produced sufficiently often in the RG for this to happen. 
Moreover, the extent to which such isolated sites can be 
associated, as in one dimension, with domain walls be- 
tween topologically distinct regions is also not clear. An- 
other obvious question is whether any of the phase transi- 
tions mentioned above are controlled by infinite-disorder 
fixed points. Although we have not studied all these ques- 
tions carefully, it is clear that our strong-disorder RG 
approach, implemented numerically, provides the natu- 
ral tool for such investigations in more than one dimen- 
sions. [So far, we have extensively studied only a related 
2d bipartite ijandom hopping problemEj using a similar 
RG approacht§ and found that such Griffiths effects do 
occur in this case: When we introduce enforced dimer- 
ization into this random hopping problem, we observe a 
transition between different Griffiths phases and tenta- 
tively conclude that the critical point is controlled by an 
infinite-randomness fised point. Details of this study will 
be reported elsewherdij and could be relevant for certain 
models for spinless T-invariant 2c? superconductors in the 
presence of disorder.] 
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